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"As  the  matter  becomes  clear,  one  tends  to  forget  the  difficulties  and  the  mystery  which  sur 
round  the  problem  in  the  beginning.  What  once  was  strange  becomes  first  evident  and  then  natural 
The  problem  ceases  to  be  a  problem;  it  fades  into  nothing  and  disappears  in  the  night  of  the  past.” 

V.  Kourganoff 
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RADIATIVE  HEAT  TRANSFER  AND  SELECTIVE  EMISSION  IN 
PLANAR  SLAB  COATINGS  ON  HOT  EMISSIVE  SUBSTRATES 


INTRODUCTION 


This  study  determines  the  theoretical  conditions  necessary  for  selective  emission  of  infrared  (IR) 
radiation  from  a  smooth  planar  coating  in  thermal  contact  with  a  hot  substrate.  The  standard  for  ther¬ 
mal  radiation  output  is  the  blackbody,  which  absorbs  and  emits  more  thermal  energy,  either  total  or  in 
an  arbitrary  spectral  interval,  than  any  other  type  of  source  at  the  same  temperature.  The  spectral 
distribution  curve  of  a  blackbody  provides  the  limiting  envelope  for  the  other  type  of  radiators.  The 
figure  of  merit  as  to  how  a  real  source  compares  to  a  blackbody  is  known  as  the  emissivity  ex  and  is 
defined  bv 


yy  source 


where  w™urce  denotes  the  spectral  emissive  flux  {W  cm  2  gm  ')  of  the  source  at  wavelength  X  and 
W*  denotes  the  spectral  blackbody  emissive  flux  that  is  given  by  the  Planck  equation  [1] 


(f)  = 


X5  eCl'XT-l  ' 


where 


h  =  (6.6256  ±  0.0005)  x  10-34  W  s 2  is  Planck’s  constant, 

T  is  absolute  temperature,  K , 

c  =  (2.997925  ±  0.000003)  x  10  10  cm  s"1  is  the  velocity  of  light, 
c,  =  2rhc2  =  (3.7415  ±  0.0003)  x  104  W  cm-2  f tm 4  is  the  first  radiation  constant, 
c2  =  ch/k  =  (1.43879  ±  0.00019)  x  104  gm K  is  the  second  radiation  constant, 
kh  =  (1.38054  ±  0.00018)  x  10~23  Ws  K~x  is  Boltzmann’s  constant,  and 
X  is  the  wavelength,  /xm. 

From  Eq.  (1),  three  classes  of  sources  can  be  distinguished  by  the  spectral  emissivity 
•  blackbody  ex  equals  1 , 


•  graybody  rx  is  a  constant  less  than  1 ,  and 

•  selective  emitter  €x,  a  variable  dependent  on  X. 
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A  quantitative  model  to  predict  the  1R  output  from  isothermal  coatings  in  contact  with  various 
substrates  can  be  formulated  from  the  solution  of  the  radiative  transfer  equation.  The  radiative 
transfer  equation  is  a  first-order  differential  equation  that  governs  the  spatial  distribution  of  the  spec¬ 
tral  radiance  in  an  absorbing  and  emitting  medium  The  radiative  transfer  equation  is  well  known  to 
the  heat-transfer  engineer  and  the  astrophysicist.  However,  in  most  engineering  applications  the  net 
flux  and  temperature  distribution  within  the  system  is  of  primary  interest.  A  common  engineering 
approximation  is  to  treat  the  medium  as  a  graybody  and  compute  the  net  radiative  flux  and  tempera¬ 
ture  distribution  under  this  assumption.  There  has  been  very  little  formal  analysis  of  spectral  emis¬ 
sion  in  such  systems. 

In  astrophysics,  it  is  recognized  that  the  absorption  coefficient  of  a  stellar  atmosphere  varies  not 
only  with  radial  position  but  also  with  wavelength.  Nevertheless,  many  of  the  theoretical  analyses  of 
astrophysical  observations  define  a  mean  absorption  coefficient  and  treat  the  stellar  atmosphere  as  a 
graybody  for  mathematical  simplicity.  Since  our  concern  is  with  spectral  emission,  the  analysis  for¬ 
mulated  here  must  avoid  the  graybody  approximation.  In  the  next  section  a  formal  solution  for  the 
spectral  IR  output  of  an  isothermal  composite  slab  is  obtained  from  the  radiative  transfer  equation. 

To  make  the  analysis  mathematically  tractable  we  introduce  a  number  of  assumptions.  The 
materials  of  the  composite  slab  are  considered  homogeneous  so  that  the  attenuation  of  radiation  by 
scattering  can  be  neglected  compared  to  absorption.  In  addition,  the  materials  are  considered  inotro¬ 
pic  so  that  the  medium  can  be  characterized  by  a  single,  uniform  index  of  refraction.  Local  thermo¬ 
dynamic  equilibrium  (LTE)  is  assumed  to  exist  throughout  the  medium  so  that  Kirchhoffs  and 
Planck's  laws  are  valid.  The  formal  solution  obtained  in  the  next  section  is  used  to  investigate  the 
effects  of  various  input  parameters  on  the  IR  output  of  the  composite  slab.  The  effects  of  coating 
thickness,  coating  absorption  coefficient,  substrate  eniissivity.  and  temperature  are  analyzed. 

ENERGY  TRANSFER  EQUATIONS 

The  general  equations  for  a  conducting  and  radiating  medium  can  be  derived  by  making  an 
energy  balance  on  an  arbitrary  volume  of  matter  This  energy  balance  equates  the  time  rate  of  change 
of  the  energy  stored  in  the  volume  element  to  the  sum  of  the  net  heat  flow  into  the  element  resulting 
from  conduction,  radiation,  and  heat  generation  within  the  volume  element.  The  resulting  equation 
for  the  temperature  distribution  is 


pc),  dT{rdt  ~  =  -  V  •  (£.  +  qr)  +  S(r ,  /), 
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is  the  density . 

is  the  specific  heat  at  constant  pressure, 
is  the  thermal  conductivity,  anil 

is  the  energy  generation  rate  per  unit  volume  per  unit  time. 
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Equation  (4)  is  a  partial  differential  equation  for  temperature  T .  Once  the  temperature  distribu¬ 
tion  is  known,  the  heat  transferred  by  conduction  is  readily  determined  from  Eq.  (5a).  Since  the  radi¬ 
ation  term  in  Eq.  (5b)  depends  not  only  on  the  local  temperature  but  also  on  the  entire  surrounding 
radiation  field,  the  energy  equation  is  an  integrodifferential  equation  for  the  temperature  distribution 
in  the  medium.  The  divergence  of  the  radiative  flux  represents  the  net  rate  of  gain  for  radiant  energy 
per  unit  volume  at  any  point  P(r)  resulting  from  the  excess  of  absorbed  over  emitted  radiation.  The 
expressions  for  the  radiative  flux  and  its  divergence  are  discussed  below. 

Radiation  Transport  Equations 

Consider  the  radiative  transfer  in  a  plane  layer  of  material,  as  shown  in  Fig.  2,  in  which  the 
temperature  depends  only  on  the  coordinate  in  the  c  direction.  The  medium  is  of  infinite  extent  in  v 
and  y  directions,  and  the  boundary  conditions  are  such  that  the  temperature  and  radiation  fields  do  not 
depend  on  x  and  y.  The  layer  is  an  absorbing,  emitting,  and  scattering  planar  slab  of  thickness  d 
that  is  in  contact  w  ith  a  hot  substrate  that  is  maintained  at  a  uniform  temperature  Ts .  For  this  planar 
system  with  uniform  boundary  conditions,  the  radiative  heat  transfer  is  one-dimensional  and  depends 
on  the  depth  only.  The  transfer  equation  (derived  in  Appendix  A)  becomes 

dNx(s,U,  t)  a{  » 

— Js -  +  «  +  «x)  Nx(s,Q,t)  =  f((s ,  T)  +  —  |fi,  =  4tP(Q'  •  0)  Nk(s,  U’)dW  .  (6) 

where  the  spectral  radiance  due  to  emission  is  f\(s,  T). 


Fig.  2  —  Physical  model  and  nomenclature 
for  a  planar  slab  coating  on  a  hot  subsirate 


In  Eq.  (6)  a“  and  ax  are  the  spectral  absorption  and  scattering  coefficients,  respectively,  in  units  of 
cm  s  is  the  length  measured  along  an  arbitrary  direction  U.  t  is  time,  8  is  the  polar  angle  as  dep¬ 
icted  in  Fig.  2,  ft  is  the  solid  angle  around  the  direction  of  propagation  ft,  and  P(ft'  ft)  is  the  phase 
function.  Equation  (6)  assumes  that  the  participating  medium  is  isotropic  and  homogeneous,  other¬ 
wise  a,"  and  ay  would  he  functions  of  direction  and  position.  Since  r  =  ,v  cos  8,  the  directional 
derivative  can  be  expressed  in  terms  of  the  derivatives  w  ith  respect  to  the  space  coordinate  :  as 
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where  p  is  the  cosine  of  the  angle  6  between  the  direction  6  of  the  oeam  and  the  ;  axis.  It  is  con¬ 
venient  to  define  the  optical  thickness  rx: 

Tx  =  lo  (c*x  +  "x)  d  ’  (8) 

The  independent  variable  s  in  Eq.  (6)  can  be  rewritten  in  terms  of  rx. 

MODEL  I:  Nonscattering  System  at  Uniform  Temperature  Ts 

For  a  homogeneous,  isotropic  coating  that  absorbs  and  emits  radiation  but  does  not  scatter  it 
(ax  =  0),  Eq.  (8)  reduces  to  tx  =  axz  and  the  radiative  transfer  equation  is  written 

aNx(rx,  p,  <p)  /x(rx,  T) 

M  - 7 -  +  Nx(tx,  p,  <fi)  =  - - - ,  (9a) 

Stx  a° 


with  the  boundary  condition  at  rx  =  0 


Nx(0,  p)  =  exNx  (Ts ) , 


where  ex  is  the  substrate  emissivity  and  the  reflected  energy  from  interfaces  at  rx  =  0  and  rx  = 
are  assumed  here  to  the  negligible.  For  azimuthal  symmetry,  the  radiance  is  independent  of  and 
Eq.  (9a)  becomes 


a/Vx(rx,  p)  ./x(rx,  T) 

M - 7 -  +  Nx(rx,  p)  =  - - - 

3rx  «xfl 


0  <  rx  <  rw  . 


The  spectral  radiance  Nx(tx,  p)  in  the  absorbing  and  emitting  coating  layer  is  determined  from  the 
solution  of  Eq.  (10).  Before  arriving  at  the  formal  solution  it  is  desirable  to  divide  Nx(rx,  p)  into  two 
contributions:  the  radiance  directed  in  the  positive  direction  (p  >  0)  denoted  by  Nf  (rx,  p)\  and  the 
radiance  directed  in  the  negative  direction  (p  <  0)  denoted  by  Nf(rx,  p).  Equation  (10)  is  then 
separated  into  Nf  (tx,  p)  and  Nf  (tx,  p)  components,  and  the  resulting  equations  along  with  the  asso¬ 
ciated  boundary  conditions  are  given  as 


^x+(rx,p)  .  .i+,  %  ft(rx,T) 

M - 7 -  +  Nx  (rx,  p)  =  - - - 

dT x  «x 


Nf  (0,  p)  =  eX(Ts), 


P  >  0, 


(Tx.  M)  x  fZ(TX,  T) 

' — 37 —  +  <Tk-  ^ - tt-  ' <  "• 


^x  (rw.M)  =  r2\H\  (Txa.  -/*). 
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where  r2x  denotes  the  specular  reflectivity  off  the  boundary  surface  7X  =  T>J.  The  reflectivity,  r2X. 
can  be  computed  for  a  smooth  surface  from  the  Fresnel  equations.  Using  an  integrating  factor  and 
then  combining  the  two  terms  on  the  left-hand  side  of  Eq.  (11a)  into  a  total  differential  yields 


r/(/Vx+  exp  (rx/fi)] 


— ~  /x(tx,  T)  exp  (tx//i)  h  >  0. 


(13a) 


In  a  similar  manner,  Eq.  (12a)  becomes 


d[Nx  exp  (rx/g)] 
d  Tx 


- -  /jf(rx,  T)  exp  (rx//i)  g  <  0. 


(13b) 


Integration  of  both  sides  of  Eqs.  (13a)  and  (13b)  with  respect  to  rx.  and  rearranging  terms  using  the 
boundary  conditions  gives 


^  .  exp(-7X/g) 

AY  <7X.  Ml  =  (0,  n)  exp  ( —  rx//i)  +  -  f  '/l  (rx.  7)  exp  (t{//z)  drx.  (14a) 

0 

and 

exp(-7X/g)  r 

AY  (TX,  m)  =  wx  (tw.  /a)  exp  [(rw  -  7X)/g]  +  - - - (  /x  (rx.  7)  exp  (rx/g)  rfrx.  (14b) 

Equations  (1 4a)  and  (14b)  are  formal  solutions  for  the  spectral  radiance  in  the  positive  and  nega¬ 
tive  ;  directions,  respectively.  In  general  the  temperature  7  is  a  function  of  rx.  However,  if  the 
coating  is  thin,  temperature  gradients  can  be  neglected  and  by  using  the  isothermal  assumption.  7 
may  be  approximated  by  7S .  With  this  approximation  j7x(rx,  7()  is  a  constant  and  can  be  factored 
from  inside  the  integral  sign  and  Eqs.  (14a)  and  (14b)  are  rewritten  as 


'fe(  7  )  r 

A V<rx,  ft)  =  AV  (0,  ft)  exp(-rx/(i)  +  — ■ -  — '  ■  f  ‘  exp  (fry  -  rx)//st)]  dr{ 


(15a) 


/vx  (TX.  /X)  =  Wx“  (rw.  M)  exp  [(7^/  -  rx)/#t|  -  --  f  'Jexp[(7i  -  rx)//x]  (15b) 

^ax  r> 

or,  by  replacing  fi  by  -g.  Eq.  (15b)  becomes 


A'x  (7X.  -M)  =  Arx“  (rw, 


-/*)  exp  [  -  (7^,  -  tx)/h]  + 


Y)/mWtx.  (16) 


Assuming  that  the  medium  is  in  local  thermodynamic  equilibrium  and  that  Kirchhoffs  law  is  valid, 
the  spectral  radiance  resulting  from  emission  is  given  by 

.K<TS)  =  (T’J-  (17) 
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To  determine  the  net  spectral  radiance  at  the  optical  thickness  7X,  Eqs.  (15a)  and  (16)  must  be  com¬ 
bined  as 

(Vf‘  (7X,  ft)  =  Nx  (TX,  ft)  -  Nx  (7X,  ft) .  (18) 

Integrating  Eqs.  (15a)  and  (16)  with  respect  to  rx  and  using  the  thermodynamic  equilibrium  condition 
from  Eq.  (17)  gives 

(t\,  h)  =  Nx+  (0,  n)  exp  (  7X/ ft)  +  N£(TS)  [1  -  exp  (-7X/g)].  (19a) 

Nx  (tx,  ~fi)  =  Nx  ( Txd ,  ~fi)  exp  [(7X  -  r^)//i]  -  n£(Ts)  [exp  j(7X  -  r^)/g)  -  1],  (19b) 

By  suing  the  boundary  conditions  for  A/x+( 0,  ft)  and  A/x_  (7^,  -ft)  from  Eqs.  (lib)  and  (12b)  the 
above  equations  become 

( Tx ,  n)  =  (Ts)  exp  (  7X/ ft)  +  (Ts)  [1  -  exp  (-7X)/g]  (20a) 

Nx  (tx<  -m)  =  r2x[((Nbx(Ts)  exp  (-7^/m)  +  N$(TS)  (1  -  exp  (-T^/ft)))  exp  [(7X  -  7^)/^] 

-  (V£(7;)[exp  ((7X  -  7^)/Mj  -  1],  (20b) 

Substitution  of  Eqs.  (20a)  and  (20b)  into  Eq.  18  gives 

Nx*(tx,  ft)  =  elNbx(Ts)  exp  (-7x//i)  +  N^(TS)  [1  -  exp  (  — 7X/^)] 

-  r2x  [txN\(Ts)  exp  ( ~tu/h )  +  N[(TS)  (1  -  exp  (-7 u/ji))]  exp  [(7X  -  TU)/ft ] 

4-  Nbx(Ts)  [exp  |(7X  -  r^)/ft)  -  1],  (21) 


The  spectral  radiance  observed  at  the  surface  of  the  slab  at  normal  incidence,  i.e.,  6  =  0°  (ft  =  1)  is 
determined  by  setting  7X  =  t xj  in  Eq.  (21)  that  then  reduces  to 


Nxe'(Txd)  =  Nbx(Ts)\  1  -  r 2X]  [exJ  exp  (-7^)  +  1  -  exp  (-7W)]. 


(22) 


Equation  (22)  is  the  central  result  of  our  analysis.  This  is  the  general  solution  for  the  radiance 
emitted  from  an  absorbing,  emitting,  nonscattering  planar  slab  coating  of  thickness  cl  on  a  hot  sub¬ 
strate  of  emissivity  fx.  The  physical  significance  of  the  various  terms  in  Eq.  (22)  are  as  follows.  The 
first  term  on  the  right-hand  side  is  the  contribution  from  the  hot  substrate  attenuated  by  the  optical 
thicknes'  the  second  term  is  the  contribution  to  the  radiance  from  the  emission  along  the  path 

f  ’  i  ■  >  -  7,  :  »he  factor  it  r,  \  account-  for  the  specular  reflect  is  its  at  the  boundan  surface 

h-.'i'A  c'i  ih  ■  o.uiiv’  and  die  outside  m,-  (  mod.  i.its’n  eh.-d’  >'n>\  be  ebon  to  s- >me  special  esc 

•  die  mpiit  parameters 
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Case  I:  Substrate  Is  a  Blackbody  (ex  =  1  and  r2x  =  0). 


Equation  (22)  reduces  to 


This  leads  to  Conclusion  I: 


=K(L)[  1]  =<(TS). 


If  a  coating  of  a  selective  emitter  is  applied  to  a  blackbody  substrate  and  forms  a  non¬ 
scattering,  absorbing,  and  emitting  layer  in  thermodynamic  equilibrium  with  the  sub¬ 
strate.  then  the  composite  system  can  only  emit  as  a  BLACKBODY. 

Case  II:  Role  of  Optical  Thickness 

The  optical  thickness  r ^  was  defined  by  Eq.  (8)  and  related  the  absorption  coefficient  ax  and 
the  geometric  thickness  d  of  the  layer.  In  an  optically  thick  regime,  the  quantity  r ^  approaches 
infinitv.  The  exponential  terms  in  Eq.  (22)  approach  zero  and  when  r2\  =  0,  this  equation  reduces  to 

A\ne,Ow)  =  N\(TS)  [1]  =  Nhx(Tx).  (24) 

This  result  leads  to  Conclusion  II: 

In  an  optically  thick  limit  the  composite  system  described  by  Fig.  2  can  emit 
only  cis  a  BLACKBODY 

This  conclusion  is  independent  of  the  nature  of  the  substrate  and  the  selective  emission  charac¬ 
teristics  of  the  coating.  It  is  also  valid  for  emission  from  hot  gases  such  as  those  found  in  IR  flares 
1 2 1 .  The  addition  of  selective  emitters  to  a  hot  flare  plume  cannot  change  the  spectral  IR  output  if  the 
plume  is  optically  thick.  Figure  3  shows  an  example  of  how'  the  optical  thickness  effects  the  IR  radi¬ 
ance  of  sheets  of  ordinary  window  glass  of  various  thicknesses  at  1000°C.  Note  that  the  spectral 
emissive  power  approaches  that  of  a  blackbody  as  the  sheet  thickness  r  goes  to  infinity  for  nonzero 
spectral  absorption  coefficients.  The  reflectivity  factor  in  Eq.  (22)  causes  the  emissivity  to  differ 
from  1.0  as  the  glass  thickness  approaches  infinity. 


When  < 


when  r2\  =  0.  Eq.  (22)  reduce; 


=  0,  the  optical  thickness  tx  is  0.  the  exponential  terms  approach  1  in  Eq.  (22),  and 
Fn  f 22 1  reduces  to 


1,'nder  these  circumstances  the  spectral  radiance  is  determined  entirely  by  the  substrate  emissivity  and 
temperature. 

The  spectral  transmissivity  I\  measured  in  the  direction  of  piopagation  is  defined  as 

Tx  =  e  T'  =  e  <:  .  (26) 
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The  conservation  of  energy  within  the  coating  is  expressed  by 

ax  +  I\  +  rx  =  1,  (27) 

where  ax  is  the  spectral  absorptivity  and  rx  is  the  spectral  reflectivity. 

In  the  present  analysis,  the  spectral  reflectivity  is  assumed  to  be  negligible:  i.e.,  rx  •»  0  and 
Kirchhoff  s  law  is  assumed  to  be  applicable,  which  implies  that  ax  =  ex.  Equation  (27)  then  can  be 
written  in  terms  of  the  emissivity  of  the  coating  layer  of  thickness  d  as 

a{  =  1  -  e~ayd  =  1  -  exp  (-r^)  =  cx.  (28) 

Substituting  Eq.  (28)  into  Eq.  (22)  yields  the  general  solution 

=  N((TS)  [e{  exp  (-rw)  +  e{],  (29) 

In  the  optically  thin  regime,  approaches  0,  and  Eq.  (29)  is  approximated  by 

Nr  (rw )=K(Ts)[e{  +  -  e{e£l ,  (30) 

where  ex  is  the  emissivity  of  the  coating.  This  equation  expresses  the  net  spectral  radiance  in  terms 
of  measurable  or  assumed  spectral  emissivities. 

Case  III:  Opaque  Substrate  of  Negligible  Emissivity  (tx  <  <  1) 

If  the  substrate  emissivity  rx  in  Eq.  (29)  is  small,  the  intensity  equation  reduces  to  the  approxi¬ 
mate  form 

Nx'W  =  Nbx(Ts)  Uxl-  (3D 

Figure  4  shows  the  measured  selective  emission  of  a  finite  thickness  of  silicon  dioxide  (tx  =  0). 
Equation  (31)  reveals  that  if  a  selective  emitter  is  applied  to  a  substrate  of  much  lower  spectral  emis¬ 
sivity  (e  g.,  a  metal),  the  radiated  intensity  from  the  composite  system  is  dominated  by  the  emissivity 
of  the  coating  layer. 

This  leads  to  Conclusion  III: 

If  a  coating  of  a  selective  emitter  is  applied  to  a  metallic  (low  emissivity,  i.e., 

IR  blocker)  substrate  and  forms  a  nonscattering,  absorbing,  emitting  layer  in 
thermodynamic  equilibrium  with  the  substrate,  the  spectral  radiance  of  the 
composite  system  will  be  dominated  by  the  EMISSIVE  PROPERTIES  OF 
THE  COATING 

As  an  example  of  this  case,  consider  a  hot  steel  substrate  at  1273  K  (1000°C)  coated  with  a 
lave:  of  silicon  dioxide  in  thermal  equilibrium.  Figure  5  shows  the  normal  spectral  emissivity  for 
these  materials  a'  the  above  temperature.  T  he  radiance  lor  a  blackbody  and  the  steel  and  silicon  diox- 
i  cmi'ooie  1 1 .  tli.  2  to  1 5-/;!!';  region  at  this  temperature  can  be  compared  by  using  hq.  (3Ui.  Fig¬ 
ure  (i  shows  the  results,  which  demonstrate  the  selective  IR  emission  lor  this  svstem. 
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fig.  4  Comparison  of  the  1R  output  of  Si0;  1 1 10  mils  thick) 
and  a  blackbodv  at  1000'C 


big  3  Normal  spectral  emissivitv  of  SiO:  1 1  10  mils  thick) 
and  steel  at  IOOO ’C  |SiO;  data  taken  from  Ref  4  and  steel 
data  taken  from  Ref  5 1 
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Appendix  A 

DERIVATION  OF  THE  RADIATIVE  TRANSFER  EQUATION 


The  propagation  of  radiation  at  any  point  in  a  medium  cannot  be  represented  by  a  single  vector 
as  in  the  case  of  heat  flow  by  conduction.  To  specify  the  radiation  incident  at  a  given  point  within 
the  medium,  one  must  know  the  radiation  from  all  directions  because  radiation  beams  are  independent 
of  one  another.  The  amount  of  radiation  energy  transmitted  by  the  ray  in  any  given  direction  per  unit 
time  is  the  spectral  radiance.  To  define  the  quantity  we  consider  an  element  of  surface  area  dA 
located  about  the  space  coordinate  7  as  illustrated  in  Fig.  Al.  Let  ti  be  a  normal  unit-direction  vector 
to  dA  and  let  dEx  denote  the  amount  of  radiative  energy  in  the  wavelength  interval  X  to  X  +  d\,  con¬ 
fined  to  an  element  of  solid  angle  d ft  around  the  direction  of  propagation  ft  streaming  through  the  ele¬ 
mental  surface  dA  during  the  time  interval  between  r  and  t  +  dt .  Let  0  be  the  polar  angle  between 
ti  and  the  direction  of  propagation  ft.  The  spectral  radiance  /Vx  is  defined  as 


dA  cos  0  dt  dU  d\ 


dfl 


Fig.  Al  —  Symbols  for  the 
definition  of  radiance 


The  spectral  radiance  in  a  given  direction  in  a  nonattenuating  and  nonemitting  medium  with  constant 
properties  is  independent  of  position  along  that  direction.  This  invariance  of  spectral  radiance  when 
no  attenuation  or  emission  is  present  provides  a  convenient  way  of  specifying  the  magnitude  of  any 
attenuation  or  emission  as  these  effects  are  given  directly  by  the  change  in  radiance  with  distance 
traveled  through  the  medium.  In  the  engineering  community,  the  quantity  defined  by  Eq.  (Al)  is 
often  called  the  spectral  intensity. 

Radiation  traveling  along  a  path  is  attenuated  by  absorption  and  scattering  and  enhanced  by 
emission  (both  spontaneous  and  induced)  and  radiation  scattered  in  the  path  direction  from  other 
direct  ions  vs  i  t  h  i  n  the  medium.  The  processes  of  absorption,  emission,  and  scattering  can  be  written  as 
an  energy  balance  over  a  small  volume  element  of  the  medium.  Consider  an  elementary  cy  linder,  the 
avis  of  which  is  directed  along  a  given  direction  of  propagation  ft  (cos  0=1).  Let  the  area  of  the 
evhnder's  base  be  dA  ,  and  let  the  cylinder's  height  be  d.\  .  The  spectral  radiance  where  the  ray  enters 
the  c\  linder  is  denoted  by  .V . .  and  the  spectral  radiance  leaving  the  vvlindcr  i>  denoted  hv  .V  *  </,V 
i  me  meoium  within  the  evlindet  interacts  with  the  radiation,  the  net  change  in  the  spectral  radiance 
dA  .  H  -  equation  ;  AI  \  the  quantity 
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<1\\  tlA  dt  dtt  d\ .  i  A2) 

represents  the  difference  in  the  radiative  energy  crossing  the  surface  JA  at  .s  -r  ds  and  s  in  the  time 
interval  dt  about  t  and  in  the  wavelength  interval  d\  about  X,  and  contained  in  an  element  of  solid 
angle  Jfi  about  the  direction  ft.  Let  Gx  represent  the  net  gain  of  radiative  energy  by  the  beam  in  the 
volume  element  dAds  per  unit  volume,  per  unit  time  about  r,  per  unit  wavelength  about  X,  and  per 
unit  solid  angle  about  the  direction  ft.  Then  the  quantity 

G^dAdtdQdXds  ,  (A3) 

represents  the  net  gain  of  radiative  energy  by  the  beam  contained  in  an  element  of  solid  angle  Jft 
about  the  direction  ft,  in  the  time  interval  dt ,  in  the  wavelength  interval  r/X,  and  in  the  cylindrical 
volume  element  dAds .  By  equating  the  quantities  (A2)  and  (A3),  we  obtain 


The  spatial  derivative  in  Eq.  (A4)  assumes  that  the  observer  is  moving  with  the  radiation  beam  at  a 
speed  c  (Lagrangian  frame).  The  distance  ds  traversed  by  the  beam  in  time  dt  is  given  by 

ds  =  cdt  .  (A5) 

To  compute  the  derivative  with  respect  to  space  in  a  fixed  coordinate  system  (Eulerian  frame),  one 
must  use  the  substantial  derivative  and  Eq.  (A4)  becomes 

|  3NX 

-  +  ft  •  V  N\  —  Gx  (A6) 


If  the  direction  of  propagation  is  along  the  path  5.  Eq.  (A6)  can  be  written  as 

1  <3jVx  dNx 

- aT  +  ~T~  =  (A7) 

<'  of  ds 

In  most  applications,  the  first  term  in  this  equation  can  be  neglected  because  of  the  large  magnitude  of 
the  speed  of  propagation  c.  We  now  derive  explicit  expressions  for  the  various  components  of  Gx: 


G  =  G“  +  G'\  +  Gh 


+  GfM 


where  Gx  is  absorption.  Gx'  is  out-scattering.  Gx  is  in-scattering  .  and  GXM  is  emission. 

Consider  radiation  of  radiance  N(s.Q)  traveling  in  an  absorbing,  emitting,  and  scattering 
medium.  As  the  radiation  passes  through  the  cylindrical  volume  element,  it  is  attenuated  bv  ahsorp 
tion.  The  decrease  in  energy  per  unit  volume,  per  unit  time  about  t .  per  unit  wavelength  about  X. 
and  per  unit  solid  angle  about  the  direction  fl  is  given  by  the  phenomenological  Bouger  Lambert  Law 

Gx  =  -  <Nk.  (Ad) 


-r.  A 


'iWJirJrSw.v 
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The  spectral  absorption  coefficient  uj'  is  a  physical  property  of  the  medium  and  has  the  units  ot 
reciprocal  length.  In  general,  it  is  a  function  of  temperature,  pressure,  material  composition,  and 
radiation  wavelength.  The  absorption  process  considered  in  fq  (Adi  includes  both  "true  absorption" 
and  stimulated  emission  because  both  processes  are  proportional  to  the  radiance  ot  the  incident  radia 
non  The  net  absorption  of  the  medium  may  be  positive  or  negative,  depending  on  whether  "true 
absorption"  or  stimulated  emission  dominates. 

The  extent  of  scattering  within  the  medium  depends  on  the  number  density  of  scattering  particles 
i.\„  i  within  the  medium  and  the  scattering  cross  section  iu“)  ot  an  individual  scattering  particle.  1  he 
cross  section  for  scattering  is  the  apparent  area  that  a  particle  presents  to  an  incident  beam  insofar  as 
the  particle's  ability  to  deflect  radiation  from  the  beam  is  concerned  The  scattering  cross  section 
depends  on  the  particle  si/e.  geometry,  material  composition,  radiation  wavelength,  and  polarization 
I  he  decrease  in  spectral  radiance  per  unit  time,  volume,  solid  angle,  and  wavelength  resulting  from 
the  scattering  out  of  the  incident  beam  direction  is  given  by 


<C  V„.\\ 


where  denotes  the  spectra!  scattering  coefficient  The  scattering  coefficient  can  be  regarded  as  the 
reciprocal  of  the  mean  free  path  that  the  radiation  traverses  before  being  scattered. 

Since  radiation  is  scattered  from  paths  adjacent  to  the  volume  element  under  consideration,  the 
scattering  process  can  also  enhance  the  spectral  radiance  along  the  distance  ds .  To  calculate  the 
enhancement  from  this  incoming  scattering,  the  directional  distribution  of  the  scattered  radiation  is 
required.  This  distribution  is  described  by  an  angularly  dependent  phase  function  Px[tt'  Sl).  Physi¬ 
cally.  the  phase  function  is  interpreted  as  the  scattered  radiance  in  a  particular  direction  divided  bv  the 
radiance  that  would  be  scattered  in  that  direction  if  the  scattering  were  isotropic.  If  the  scattering 
particles  in  the  medium  are  modeled  as  isotropic  spheres  and  the  medium  has  no  preferential  direction 
for  scattering,  then  the  phase  function  depends  only  on  the  angle  between  the  incident  direction  ft  and 
the  scattered  direction  O'.  Mathematically,  the  phase  function  can  be  related  to  a  probability  density. 
The  quantity  defined  by 


PAW  0)  dQ'  . 


(All) 


represents  the  probability  that  an  incident  beam  of  radiation  traveling  in  direction  0  will  be  scattered 
into  an  element  of  solid  angle  dW  about  the  direction  0'  The  scattering  of  the  incident  radiation. 
,VA(f .  f2 ) .  by  the  medium  per  unit  time,  per  unit  volume,  per  unit  wavelength  into  an  element  of  solid 
angle  t/0'  about  the  direction  12'  is  given  by 


<rx'.VA(y .  12)  </12  ~-A\(12'  12)  d  12' 

47T 


( A 1 2 ) 


In  general,  the  scattered  radiation  will  be  incident  on  the  volume  element  from  all  directions 
The  integration  of  Fq.  (A  12)  over  all  solid  angles  of  incidence  yields 


f;/;  =  ~  j/V$2  12)  ,Vxl.i .  12)  di 2. 


<  A  l  A ) 


which  represents  the  scattering  of  radiation  incident  on  the  volume  element  from  all  directions  into  the 
solid  angle  dW  about  the  direction  12'  per  unit  time,  per  unit  volume,  and  per  unit  wavelength 
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If  the  medium  is  in  local  thermodynamic  equilibrium,  spontaneous  emission  along  path  length  ds 
enhances  the  radiance  in  the  s  direction  by 

G™  =  aX  •  (A  14) 

where  N x  is  Planck's  function.  Equation  (A14)  is  a  mathematical  statement  of  Kirchhoffs  law.  The 
assumption  of  local  thermodynamic  equilibrium  simply  means  that  the  matter  contained  in  the  cylin¬ 
drical  volume  is  in  thermodynamic  equilibrium  with  itself  but  not  necessarily  with  the  radiation  field 
I A 1  ] . 

Given  these  explicit  expressions  for  the  components  of  Gx,  the  radiative  transfer  equation 
immediately  follows  from  Eqs.  (A4)  and  ( A 8 ) : 

dNx  a*  _  . 

—  =  -  K  +  axu)  +  —  (  Px  (ft'  •  ft)  Nx(s.  ft)  dft  +  c*X  .  (  A  15) 

as  4tt  J 


or 


dNx(s,Q,t)  ax  .  „  . 

- - - +K  +  ftxi)iVx(i,fi.r)  =  aX(s.n  +  —  L  Px(ft'ft)JVx(*.ftVft'.  (A16) 

as  4tt  u  =  4* 

For  further  details  see  Ref.  A2. 
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Appendix  B 

AVERAGE  (HEMISPHERICAL)  SPECTRAL  RADIANCE  CALCULATION 


The  average  spectral  radiance  for  azimuthally  symmetric  radiation  is  defined  as 


Gx(rx) 


r  r 

J0  J  ,  Wx(rx,  h)  dn  dtp 


C  \[{dVd* 


Jo 


(Bl) 


which  in  the  present  case  of  radiances  in  the  forward  and  backward  directions  becomes 


GX<  rx» 


Ox.  n)  dn 

r  I 

f  dfl 

Jl) 


0 

j  _ ,  'V  I f\  •  /*>  d n 


(B2) 


Substitution  of  Hq.  ( 15a).  and  ( 16)  into  Eq.  (B2)  yields  the  formal  solution  for  the  spectral  radiance  as 
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(i ( rv )  =  j  ,VX  (0.  n)  exp  (  ~Tx/fi )  d^ 
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By  usine  the  boundary  conditions  lor  \.  (U./,i.  anil  A  .  )  from  Eqs  ( I  I  hi  anil  <  1 2b)  and 

the  the!  mody  namk  equilibrium  condition  Irom  l  q  (  I  7 1  the  above  equation  becomes 
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LAFEMINA  AND  LADOUCEUR 


By  using  the  exponential  integral  function  defined  as 


Eq.  (B4)  is  written  as 


En(z)  =  |o  rin  2e  :/v  dij. 


Gx(tx)  =  t{Nbx(Ts)E2(TX)  +  Ni  i  (rx  -  ri)dT' 


r2x^x  (T’i  )[c\f 2<^x)  +  ^:(Txa  -  rx)  +  £:(2^  “  ^x)! 


+  ^xl7'.)!/  -  Tx>  dr'. 

The  two  integral  terms  in  Eq.  (B6)  can  be  reduced  to  a  simpler  form  by  using  an  identity  of  the 
exponential  integral  as 


N\<T,)  j  (  £|(tx  -  TX)dr{  +  E|(tx  -  rx)r/rx 


=  N\(TS  )(2  -  E2(tx)  -  E:(rw  -  tk)] 


Substitution  of  Eq  (B7)  into  Eq.  (B6)  gives  the  average  spectral  radiance  as 


Gx(tx)  =  N^iT,)  |e££2(TX>  +  2  -  E2(tx) 


E2(tm  -  tx)  -  r:x|fx  E2(tx) 


+  EdT^t  -  Tx)  ~  E2(2tu  -  TX)]{ 


The  spectral  radiance  observed  at  the  surface  of  the  slab  is  determined  by  setting  tx  =  in 
Eq  (BSi,  vshich  then  becomes 

=  A/C ( 7", )  |^Edru)  +  1  -  E2(TXil)  -  r2x|fxE2(ru)  +  1  -  E:(tu)||  (B9a) 


OVr^l  =  ,Vj(7;>  |l  -  r:x]  |1  +  E2(tw)  Ux  -1)).  (B9b) 

Under  these  circumstances  the  spectral  radiance  is  determined  entirely  by  the  substrate  enns- 
sivity  and  temperature.  When  tx  is  small,  the  exponential  function  E2{tx)  can  be  expanded  into  a 
series  expansion  approximation  as 
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Ez(t\)  =  1  -  7\  +  7 7\  +  0(tx2)  , 


(BIO) 


where  ->  =  0.577216  (Euler's  constant),  and  O(t^)  are  higher  order  terms. 

Substitution  of  Eq.  (BIO)  into  Eq.  (B9)  yields 

Gx(rw)  =  Nhx(Ts)  [1  -  r2X]  [1  +  (1  -  rx  +  7rx)  Uxs  -  Dl.  (Bll) 


which  can  be  reduced  to 


G\(tu)  =  N'X(TS »  1 1  —  r:x]  Ux(  1  -  7 ’tx)  +  7'rxl ,  (B12) 


where  -> '  =  1 

By  using  rs  =  »  ;  from  Eq.  (29).  the  average  spectral  radiance  in  the  above  equation  becomes 


Cx(tw>  =  .Vx (Tj  1 1  -  r:x]  |ex(l  -  y'({)  +  >  ,  (B13) 


or 

cX(tu)  =  Atom  -  r;x]  |fj[  +  7;  -  fX 7X] . 


(B14) 


w  here  <  i  ’<  v . 

l.’pon  comparison  ot  net  radiance  normal  to  the  surface  of  the  slab  (Eq.  30)  and  the  average 
radiance  o\er  the  hemisphere  (Eq.  BI4i  with  r:K  =  0.  it  is  seen  that  the  effect  of  averaging  comes  in 
by  replacing  the  emissivity  of  the  coating  layer  by  hemispherical  emissivity  . 


